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In Ising model on the simple cubic lattice, we describe the inverse temperature j3 in terms of the 
bare-mass M and study its critical behavior by the use of delta expansion from high temperature or 
large M side. In the vicinity of critical temperature /3 C , the expansion of /3 in M has j3 c as the first 
term and M~ 1 ^ 2l/ as the leading correction. The estimation of f3 c in 1/M expansion is confronted 
with the leading and higher order corrections, even delta expansion is applied and the critical region 
emerges. To improve the estimation status of f3 c , we try to suppress the corrections by adding 
derivatives of /3(M) with free adjustable parameters. By optimizing the parameters with the help 
of the principle of minimum sensitivity which are maximally imposed in accord with the number of 
parameters, estimation of j3 c is carried out and the result is found to be in good agreement with the 
present world average. In the same time, the critical exponent v is also estimated. 

PACS numbers: 11.15.Me, 11.15.Tk, 64.60.Bd, 64.60.De 



I. INTRODUCTION 



Since the invention of lattice field theories, the border between condensed matter models and field theoretic models 
r— n is lost and the techniques in the statistical physics have been frequently used in the field theory analysis [TJ [2]. 
Traditionary, the both systems are described in terms of j3 which indicates inverse temperature or inverse bare 
coupling constant. In the field theory side, we however notice that the lattice spacing a or the equivalent bare-mass 
O | can play the role of the basic parameter describing the models. This feature naturally appears in the large N limit 
of field theoretic models. Also for finite-TV case of non-linear sigma models, it was shown that the scaling behavior 
7^ i of the inverse bare coupling has been captured from large bare- mass expansion [3] . Then, it is a small and natural 
step to take the reverse point of view into the condensed matter models. The motivation of this work is to investigate 
y—< whether such reverse approach to the critical phenomena is effective or not. 

In the present paper, we concern with the 3-dimensional Ising model on the simple cubic lattice as a theoretical 
laboratory. Let us sketch our strategy below: As the temperature approaches to the critical one from the high 
' temperature side, the correlation length £ diverges as 

£ ->/+(l-/3//3 c )-"{l + const x (l-/3/7? c ) 8 + •••}, p < B , (1.1) 

m 

where j3 c , v, /+ and stand for the critical temperature, the exponent associated with £, the amplitude in the high 
temperature phase and the exponent of confluent singularity [3], respectively. The scaling law (1.1 1 can be rewritten 



in terms of the bare-mass M, where M is defined by the magnetic susceptibility x an d the second moment fi as 

M = — , D : space dimension. (1-2) 

A 4 



First we invert ( 1.1 1 and obtain /3 — /3 C (1 — fV v ^ » + ■ ■ •). Then, since £ ~ M l / 2 in the critical region, we have 

P = Pc- MT Mh +••', (1-3) 
where the " • • •" represents the higher order corrections. Thus the critical temperature is given by the limit 

I3 C - lim (3(M), (1.4) 



and (2u) 1 in ( 1.3[ ) is interpreted as the exponent of the leading correction. The description by M is not restricted 
to f3 but is applicable to \ 5 , specific heat and maybe others. For instance, we can express x by M via the relation 
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with /3 such that x ~^ const x M^ 1 ^ 2 ^ as M — > 0. Note that we need no information of f3 c to study the small M 
behavior of %(Af) and the estimation of 7 is unbiased. 

In the study of the scaling behavior of j3 } we confine with the high temperature phase and use the 1/M expansion. 
Apparently, we need some improvement to allow the use of the series near the critical point. As a key technique, 
we apply the so-called delta expansion method on the lattice Actually in [5], the Ising model on the square 
lattice was revisited to examine the power of the method. Though a remarkable improvement on the behavior of 
high temperature expansion was shown, the discussion there ended at the semi-quantitative level and the estimation 
task of fi c and v was not attempted. The difficulty for the accurate quantitative study of critical quantities comes 



from the corrections to the asymptotic scaling, mainly from the second term in (1.3) and even from higher orders. In 
the present paper, we tackle the problem by introducing freely adjustable parameters into the naive thermodynamic 
quantities. For /3, we instead consider 

^({ Pk } 1 x)=[Y J Pk(x^) k ]p(x), x = M- 1 . (1.5) 

fe=0 

Here, po = 1 and p\ = p, p^ = a, ■ ■ ■ are adjustable free parameters. By exploiting a special property of the delta 
expanded tp, we make approximate cancellation of corrections to /3 C by seeking optimal values of parameters order by 
order. On the basis of the high temperature or large bare-mass expansion, we then try to compute /3 C and v. The 
values of /3 C and v have been computed by various methods including Monte Carlo simulation, field theoretic methods 
centered on the renormalization group and series expansions (See, for a review [7] and recent researches [BHT3]). The 
best estimation of /? c known to our knowledge is f3 c = 0.2216546(10) [T5] and f3 c — 0.2216595(15) [16] . In this work, 
we however quote modest one 

j3 c = 0.22165, (1.6) 

to which all recent literatures agree up to the last digit. For v we refer 

v = 0.630, (1.7) 

which is also agreed by all recent works within 10~ 3 order. 

This paper is organized as follows: In the next section, we briefly review the method of delta expansion on the 
lattice. In the third section, we attempt to estimate critical quantities, f3 c and z/, of the Ising model on the simple 
cubic lattice by the use of the delta expansion on Studies on the magnetic susceptibility and specific heat, including 
also for the square lattice model as a bench mark of the method, are now under the progress. To avoid a report too 
long, we hope to carry those studies forward to another publication. In the present paper, we like to present the 
essence of the idea and technical details of our approach by focusing on the relation of (3 and M on the cubic lattice. 
The conclusion is stated in the last section. 



II. DELTA EXPANSION 



To access the critical region, we make an attempt to dilate the region itself around M = 0. Given a thermodynamic 
quantity /(M), we consider the dilated function /, f(M,S) = f(M(l — 8)) with < 8 < 1. Setting the value of 8 
close to 1, the critical region, the neighborhood of M = 0, is enlarged to the region far from the origin. Then, if 
the large M series of /(M(l — 8)) still effective in some region of M can be available in the 8 — > 1 limit, the series 
may recover the original critical behavior within there. To obtain such an effective large M series of /(M, 8), we have 
found that a protocol of obtaining best large M series for f(M, 8) is to treat M _1 and 8 on equal footing: Let M in 

f(M,S) denote t~ x to avoid notational confusion. Then consider the truncated series /at(M) = X^=o a nM~ n an d 
its dilation, /iv(t _1 (l - 8)) = J2n=o a n{t/{l - 8)} n . If the full order of the scries is N, the term {i/(l - 8)} n should 
be expanded in 8 and truncated such that the sum of orders i and j is equal to or less than N where i and j denote 
respectively the orders of t and 8 [6], [17] . In this rule, {t/(l — 8)} n should be expanded in 8 up to the order S N ~ n , 
and we find 

Ml - 8)V - t» ( 1 + n8 + n ±±^8 2 + ■■■+ ^ + l)---(n + N-n-l) sN , 
V 2! (TV — n)\ J 

After the expansion in 5, we can take the limit 8 — > 1. The result gives the transform, 



{t/(l-8)} n ^C(N,n)t n , 



(2.2) 
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where 

TV! r(TV+l) 
n\(N — ny. I (n + 1)1 (TV — n + 1) 

The above transformation rule is very simple. If one has the truncated 1/M series of f(M) to order TV, one obtains 
readily the corresponding delta-expanded series. That is, given a truncated series, 

N 

f N = Y / a n M- n , (2.4) 



n=0 



the corresponding delta-expanded series -D[/jv] reads 



N 

D[f N ] = a nC{N, n)t n = f N (t). (2.5) 

n=0 

We notice that C(TV, 0) = 1 and the constant term is left invariant. 

To summarize, the delta expansion creates a new function /jv(i) associated with /jy(M) by the transform of the 
coefficient from a n to a n C(N, n) which depends on the order TV. The symbol D denotes the transformation of Jn{M) 
to /at(£). Since C(TV, n) — > ^y- (n:fixed, TV — > oo), one might think that the series becomes ill-defined in the TV — >• oo 
limit. However, in some physical models and mathematical examples, we found the evidence that within some region, 
the limit indeed exists. Consider, for example, the function f(M) = (1 + M)^ 1 = 1)™ ~ 1 M~ n . We find 

/iv(i) = X^il - !)" -1 Tt!(jv-«)! ^" = 1 ~ (1 ~ t) N an< ^ the resulting function converges to 1 for |1 — 1\ < 1 and diverges 
for |1 — t\ > 1. The point is, within the region |1 — t\ < 1, limjv-->oo /jv(i) — > 1 = /(A/ = 0). Also from general point 
of view, the result is reasonable, since formally /(Af, 1) = f(M = 0). In any way, in the present work, we assume 
that, over some region of t, the limit of the function sequence {//v(i)} tends to a constant one and 

lim f N (t)=f(M = 0). (2.6) 

iv— >oo 

In the presence of phase transition, we must deal with the case where the expansion of / around M = is not 
regular. Then we consider how in such a case the delta expansion affects the small M behavior of /(M), supposing 

that / behaves at small enough Mas/~ /(0) + fiM ai + f 2 M a2 H where < a x < a 2 < ■ ■ ■■ When t -1 (l - S) is 

substituted into M a and ((1 — 5)/t) a is expanded in 5, giving M a = T~ Q (1 — aS+- ■ •), a reasonable truncation protocol 
for the best matching with fN{t) is not found on logical grounds. Here also, we proceed along with experiences. In 



some physical models, we found that the formal extension of (2.2) 



M°>^C(N,-a)t- a , C(N,-a) = - f^* 1 ] — v , (2.7) 

1 {—a + 1)1 (TV + a + 1) 

provides us the best matching. The factor C(TV, —a) vanishes when a = 1, 2, 3, • • • and for positive non-integer case 
goes to zero as 

C(TV, -a) ^ TV- Q /T(-a + 1) ^ 0, (TV->oo). (2.8) 

Hence, the term of positive integer power of M vanishes, 

D[M a }=0, (a = l,2,3,---), (2.9) 

and the terms of fractional positive power of M decreases with the order and disappears in the TV — > oo limit. Thus, 
we define the result of delta expansion for M a by \2.7l and 



D[f] = /(0) + /iC(TV,~ ai )^ ai +..., (i»l). (2.10) 



The above two results (2.8) and rt2.9k show the main advantages of the delta expansion. From these we understand 



that the approach to the critical behavior is quicker in the 5-expanded function than the original function. 



In the present work, our task is to estimate /(0) and a.\ from the known series (2.5). In the process, we use 
derivatives of /jv- Then we remark that 

d 



t^D[f N ] = D x — f N , (2.11) 



d_ 

which states that D-operation and differentiation is commutable. It is convenient to use the following abbreviate 
notation, 



[ t Jt) klN = I ^' t^ 1 ' 2 ' 3 '-")- ( 2 - 12 ) 
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III. ESTIMATION OF /3 C AND v 



A. Preliminary study 



The Ising model on the simple cubic lattice is defined by the action 



<i,j> 



8,8 



(3.1) 



where the spin sum is over all nearest neighbour pair on the periodic lattice. Our approach is based upon the high 
temperature expansion. The magnetic susceptibility and the second moment have been computed up to (3 25 by Butera 
and Comi [TH]. From (1.2) and the result reported in [TH], we have 



/3 = 



1 6 124 312 12596 21432 
.\/ M- 3M3 ~~ M4 + 5Af 5 

148384348 797787336 17341288504 
+ 9M 9 5M 10 + 11M 11 



1330848 1745344 

7M 7 M s 
15857888272 

M 12 



2106367479672 11748802870160 263968267347944 

13M 13 " 7M 14 " + " 15M 15 
186504592354608 33924951987330804 21535692193295224 



Af 16 

4449606807205690200 

19M 19 
197756701920466780928 

7M 21 

80156432259652309452520 
23M 23 

10946582972904015563857296 



17M 17 
12821205881021198992 



M 18 



5A/20 

3442869826889278353376 
11M 22 

116948936021276297965072 



3M 24 



25M 25 



+ 0(M~ Zb ) 



£ 

n=l 



(3.2) 



The result of the delta expansion to the order N is readily obtained by multiplying nth order coefficient by C(N, n), 
giving 



JV 



P N :=D\p N ] = Y J KC{N,n)f n 



(3.3) 



n=l 



The effects of delta expansion are clearly shown in the plots of relevant functions. We show in Figure 1 the plots of 

#25) #25) #25 j #25') #2? an d #25^ I R tne nrs t g ra P n i it is implied that /?25 approaches to the correct j3 c (This is not 
the case for square Ising model). Also in the derivatives, impressive point is demonstrated: Beyond the peak about 

t ~ 0.01, #25^ shows the monotonic decreasing trend to t ~ 0.15. We have numerically checked by using the values, 



/+ ~ 0.5 for the amplitude [TH] in (1.3) and v — 0.63 (see(1.7l), the rough agreement of the behaviors between 



in i-series and its critical behavior - f -^ r C(2h,-l/2v)t- 1 l 2u (C(25, -l/2v) = 0.0170152). Same thing applies 
to the second derivatives. In addition, in the case of square Ising model, we ascertained by using known results of v 
and the amplitude that the transformed derivatives exhibit critical behaviors. Thus, we come to conclude that three 
plots in Figure 1 afford evidences emphasized in [5] that the delta expansion dilates the scaling region to the region 
of t far from the infinity. 

In the study of phase transition via series expansion technique, we rely upon heuristic assumption of the power law 
or logarithmic behaviors near transition point. Let us start the argument by supposing the scaling behavior of inverse 
temperature in general form, 



-pi 



A 2 x 



-P2 



A 3 x 



-P3 



where < pi < p 2 < P3 < ■ ■ ■ and 

P = P C - A ± C(N, - Pl )t^ - A 2 C(N, - P2 )r^ - 0(t"»). 



(3.4) 



(3.5) 
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FIG. 1: Plots of /32s, P25 , P^B > ^25' an d Z^s*- P25 ^ s shown in the gray plot and its delta expanded version p!£> in the 
black plot. In the first plot, the dotted line indicates the critical temperature f3 c — 0.22165 • ■ • which value is widely confirmed. 
In P25, the high temperature region is pushed away to the neighbourhood of the origin t = and the scaling region appears to 
cover up to point t ~ 0.15. In the second plot, we see, though ^25 is monotonic, /J^' has a peak around t ~ 0.01. This peak 

rt( 2) 

indicates the turning point from high temperature region to the scaling region. The same reasoning applies to f3 2 ^ in the third 
plot. 



The comparison of (3.4) with (1.3) gives 



Pi 



1 

2^ 



(3.6) 



and the exponents of confluent singular terms [4] enter into pi where i > 2. For instance, p% = (1 + 9)/2v. Our 
assumption is that over some region of t, 



lira ft 7v 

_/V->oo 



(3.7) 



Though the definition of /3 C is given by (1.4 1, it cannot be naively used for its estimation, because our information 
is limited to the truncated version (3.2) (One might think that Pa.de approximants of (3.2) may be useful for the 
purpose. Investigation on this direction is out of this work). Moreover, though the behavior of is remarkably 
improved as shown in the first plot in Figure 1, it is not enough yet to yield accurate estimation of /3 C . The reason 
is that the leading correction to (3 C given as AiC(N, —pi)t~ Pl is still active in the effective region of j3jy for N < 25. 
Now then, our idea upon the estimation of /3 C is to cancel dominant corrections to j3 c in (3.5 1 by introducing following 
i\) with a free parameter p, 



tp(p, x) = (3 + px 



dp 
dx 



At large x (small M), ip(p,x) reads 



(3.8) 



ip(p, x) = /3 e - y^(l - pp n )A n x p " 



n=l 



Note that for any value of p, 



lim ip(p, x) = /3 C 



Then, we deal with (5-expanded function 



4> N (p,t) = D[iP N (p,x)] = p N +pP {1) 



N ' 



where i/jn(p,x) denotes the truncated series to the order x . It follows from (3.7) that 

lim ip N (p,t) = p c . 



(3.9) 

(3.10) 

(3.11) 
(3.12) 



This should hold at any t in the supposed convergent region. We notice that, as the extension, one can build ip with 
multi-parameters. Making use of the property (3.12), we carry out estimation of /3 C and v in the following subsections. 
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B. p c and v in one-parameter case 



When dealing with i/j, the critical issue is how we should choose p value to cancel the correction in (3.5 1. The 
starting point is (3.12) which states that V>jv(p,i) approaches as N — > oo to a constant function over the conve rgent 
First, at finite N, we employ the principle of minimum sensitivity due to Stevenson [20] . From (3.12), the 

Thus, we postulate that, at a given p, the value of ipN{p,t) at the possible 
The resulting /3 C depends on p. To extract best estimation, we notice that 



region 

use of the principle is quite natural, 
stationary point gives an estimation of /3, 
if the suppression of the correction is achieved, the coefficient of the first order correction in (3.9) may almost vanish, 
leading p ~ p^ 1 — 2v. And then the stationary property of t/Sjv becomes maximal. Thus, taking reverse point of view, 
we search p which realizes the maximal stationarity and minimization of the second derivative. 

The explicit procedure goes as follows: At a given p, the estimation of f3 c is achieved at t satisfying 



5(1) 
Pn 



5(2) 
3 JV 



0. 



(3.13) 



Then, to realize maximal stationarity, we furthe r imp ose that optimal p should locally minimize the absolute va lue of 



-(21 

the second derivative \tp N | at t which satisfies (|3.13|) 
to become also the solution of the equation, 



Then, our task is to try to adjust p for the solution t of (3.13) 



^\p,t)=ffi' +pffi> = 0, 



- fl(2) 



(3) _ 



(3.14) 



— (2) 

or to make \ip N 



locally minimum. The two conditions, (3.13) and being local minimum of \i[) N |, determine optimal 
p = p" and t = t* at which the stationarity is maximally realized. In general, for a given p, the first condition (3.13) 
gives non-unique solutions of t and it is convenient to express p as a function of t, p = p(t) 
that 



From (3.13), it follows 



and then 



5(1) 
J N 
5(2)- 
3 JV 



^W) = ((4 2) ) 



4% 3) )/4 2) - 



With obtained t* and p* = p{t*), we estimate j3 c by 

/3 c (N)=^ N (p*,t* 



(3.15) 



(3.16) 



(3.17) 



We carry out this procedure from 4th order to 25th order. The order from which the real feature of our approach 
begins to show is 13th order. Here we mean by "real feature" the feature that the characteristic scaling behavio r that 



"(2) - , 

the second derivative if) N should exhibit, which we describe below: The parameter p given as a function of t as (3.15) 




0.2 0.25 



FIG. 2: p(t), ipP(p(t),t) = ((/Jjv') 2 - ^nPn^/Pn* and ipN(p*,t) at N = 13. The singular point in the upper two graphs 
comes as the solution of /3 (2) = 0. They are at t = 0.01526 and 0.14768. Note that between the two points, p is positive definite. 
The range of optimal stationary points are thus limited within 0.01526 < t* < 0.14768. We notice that the largest solution 
of = exists at t = 0.18924. As is seen in the p(t) plot, this however leads to negative p and negative p 1 1 . Hence, we 
neglect it. Two candidates of physical interest are shown by the points indicated by black triangle and gray box. Among them, 
important one is at larger t because it lies in the place where the critical region begins to appear. The lower plot shows two 
^>i3(p*,t). The gray one is for p* = 0.863705 and the black one is for p* = 0.945108. 
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should behave at t in the critical region as 

EZiPkC(N,-p k )A k t-P" 



9 



1 

Pi 



p f^~ p i A . 2 (pi-p^- {p2 - pi) +--- 



P IC(N, - Pl )A 



~(2) 

Then ip N loses t~ Pl term and has to behave as 



(3.18) 



- -(Pi -P2) 2 C(N- P2 )A 2 t-^ + 



P2 



Pi 



(3.19) 



Note that the correction vanishes in the N — > oo limit, since C(N, —p<i) —> l/r(l — p 2 ) x N~ P2 as N — > oo. At finite 
order the correction ~ const x t~ P2 remains and, therefore, we understand that the weak peak about t ~ 0.09 at 
N = 13 shows the transition from high temperature to scaling regions. This is supported by the higher order plots 

— (2) 

of ip N ( see the right graphs in Figure 3). The dotted vertical lines in the left and right graphs in Figure 2 represent 
the value of t at which /3{ 3 ' = and p — > ±oo. 

As shown in Figure 2, at 13th order, we obtain two candidates of the estimated j3 c . One comes from the zero of 

— (2) ~(2) ~ (2) 

ip\ 3 and the other from the local minimum of \ but ^13 7^ 0. The first candidate lies in the pre-scaling region 
and the other in the beginning of the scaling region. Hence, we regard the candidate at lager t as important. Actual 
values computed are /3 C (13) = 0.218638 (at t* = 0.068789) and /3 C (13) = 0.220024 (at t* = 0.116908). As we guessed, 

~(2) 

we have thus confirmed that the best choice is the one at larger t* . We note that at the point where tp N = 0, p(t) 
becomes extremal. At general order N, this is found by the direct differentiation of p(t) by t. 
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FIG. 3: p(t) (left) and ipP(t) = ((P^) 2 - P^Pn^IP^ ( ri S ht ) at TV = 24 (upper) and N = 25 (lower). The zero point of 

—(2) 

4> N (t) at smallest t > moved toward left compared to the 13t h ord er. Thus, both scaling regions have developed broader 

—(2) 

compared to the 13th order. This agrees with the assumption (3.121. At 25th order, largest solution of ip^ = exists at 
t — 0.178826. This however leads to negative p and negative pf" 1 - Hence, we neglect it. 



Then, let us turn to the high order cases, 24th and 25th orders. From the plots shown in Figure 3, we find longer 
and clearer scaling behavior than 13th order. At 24th order, it is sufficient to confine ourselves with the region p > 0. 

~(2) 

Though there are two zeroes of V'24 1 we m ay select the solution at larger t. At 25th order, the feature is similar with 
that at 13th order. The difference, however, consists in the clearness of the scaling region after the weak peak. This is 

~(2) 

because zero of "025 has moved toward the origin and the scaling region has developed. Due to the same reason with 
13th order case, we should rely upon the estimation at larger t. In this manner, we can identify the good sequence of 
estimated j3 c . The results from 20th to 25th orders are shown in Table 1. In Figure 4, the sequence of f3 c estimated 
at smaller t is shown in black triangles and that at larger t, consisting good sequence, by gray boxes. 

As the order increases, the optimized value of p gradually increases but not to be seen to converging the established 
value, 1/pi — 2v — 1.26. On the other hand, the accuracy of estimated f3 c is good for N = 25 as f3 c — 0.221117- • • 
and the relative error is about 0.24 %. The non-accurate results for p\ may be explained as follows: In our method, 



TABLE I: Estimation of f3 c and pi = l/(2v) = 0.79365 (1/pi ~ 1.26) with one parameter. 



order 


20 


21 22 


23 


24 25 




0.220933 


0.220944 0.221035 


0.221040 


0.221114 0.221117 


1/pi 


1.028607 


1.029841 1.041745 


1.042425 


1.053098 1.053408 


t* 


0.133502 


0.126393 0.135082 


0.128043 


0.136505 0.129586 




2v 



10 15 
order N 



FIG. 4: Estimated result of fi c and p 1 1 = 2u from 4th to 25th orders. Black triangle indicates the result at smaller t and the 
gray box the result at larger t. The channel of good sequence opens at 8th order for N and 13th order for odd N. 



the parameter p is selected for the total correctio ns t o be suppressed. Here notice that the dominant part of total 



corrections are composed by first few terms in (3.5) and not only by C(N, — p\)A\t~ Vl . Hence, even when the 
cancelation is successful, it does not necessarily mean that the sole information of p\ is accurately absorbed in p* , 
unless the order of 1/M expansion is extremely large. 

The sequence of good estimates starting with 8th order for even N and 13th for odd N has enough number of 
terms to do extrapolation to the infinite order by the fitting of obtained estimation. Assuming the simplest form, 
flc(N) ~ p c ~ b x N~ u , let us use the last three estimations of odd orders, 21st, 23rd and 25th. Then, we obtain 

f3 c = 0.221642, (3.20) 

with b = 0.101171 and u = 1.634371. If we use the 20th, 22nd and 24th orders results, we obtain fi c = 0.221624 with 



b = 0.10283 and u — 1.669628. Both values come to close to the established value, j3 c = 0.22165 (see (1.6)). The 
relative error is only 0.00035 ~ 0.01 %. Also for p^ 1 we have done the fitting and obtained with 21st, 23rd and 25th 
order results, 

p{ 1 w 1.2828. (3.21) 
And with 20th, 22nd and 24th orders, we have p^ 1 w 1.2640. The estimation is much improved as having relative 



discrepancy to ( 1.7[ > about 0.3 ~ 1.8 %. We can say that our extrapolation provides consistent values with the present 
world standard. 



C. /3 C and v in two-parameters case 



In this subsection, we extend our method by introducing two parameters p and a associated with fy$) and fi( 2 \ 
respectively. We expect that the problem of the inaccuracy of estimated p\ may be partially resolved by introducing 
multi-parameters. 

Consider the function i/j(p, a, x) which has the same limit as x — > oo with /?, 



ip(p,a,x) 



1 + px 



dx 



alx ioT. 



The delta expansion of ipN reads 

i) N {p,a,t) 



d \ 2 



df 



p N (t) = j3 N + p/3%> + aP 



N 



(3.22) 



(3.23) 
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Under some domain of (p, a) and some region of t, we postulate 



lim tp N (p,a,t) = P c 

N-yoo 



(3.24) 



Then, since ipN behaves at large t as, 



$ N (p,a,t) =p c - V(l - PPk +apl)C(N,-p k )A k t- p ", 



fc=i 



(3.25) 



we infer that, at least for large enough N, optimal p and a which realize the maximal stationarity may lead that the 
first two corrections in (3.25) effectively disappears. It then follows that 



1 - PPi + <JPi ~ 0, 1 — PP2 + crp 2 ~ 0. 

We use these upon the estimation of pi from optimal p and a. 

Since we have two free parameters, we impose stationarity criteria up to the second derivative, 



In terms of /?, the above condition reads as 



B (1) 
Pn 

5(2) 

Pn 



p4 2) 



^\p,a,t)=Q. 



vffi = 0, 
= 0. 



(3.26) 



(3.27) 



(3.28) 



For a given set of (p,cr), we search for the stationary points at which (3.28) is satisfied. As in the one-parameter case, 

estimated /3 C varies with values of (p, a) and we have to extract the optimal set among them. The condition is that, 

J L ~c^) 

just at t which is the solution of (3.28) , \ip N | becomes local minimum. To proceed further, it is convenient to express 
p and a in terms of the solution t for (3.28). It is then readily obtained that 



S(3)5(3) 

Pn Pn 



5(1)5(3) 

Pn Pn 



Pn Pn 



(4 2) ) 2 



A 



where 

The two parameters behave at large t as 

P 



5(2)5(4) 

Pn Pn 



(^ 3) ) 2 



(3.29) 
(3.30) 

(3.31) 



Pi +P2 

PlP2 
1 



PlP2 



+ const x r P3+P2 + 
+ const x r P3+P2 H 



(3.32) 



This shows the critical behaviors of two parameters. The critical behavior of ^ 3 \t) := ffl 3 ' (p(t) , cr(t), t) is then given 
by "0 < - 3 - > (t) ~ const x t~ P3 . The coefficient of t~ P3 denoted by "const" depends on the order N and tends to zero in the 
N — > oo limit. Therefore, as in the one-parameter case, we consider the minimization of |V^ 3 '(£)| including the case, 



(3.33) 



(3.34) 



4 3) W = 4 3) + pWf + °m { N ] = o. 

Then, with obtained t* and p* = p(t*) and a* = o~(t*), the estimation of (3 C is directly given by 

8 c (N)=i> N (p*,<r*,t*)- 

Figure 5 shows the plot of ^25 as function of t in different scales. The vertical dotted lines indicate the singularities 

of 1P25 ■ AH °f them come from zero points of A, the solutions of A = 0. The stationary point with \1p25 \ locally 
minimum but non-zero lives in the left region separated by one of the lines. In fact, the corresponding value of a 
vanishes and the solution corresponds to the one obtained in one-parameter case with smaller t* . Thus, it still lies in 
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FIG. 5: (t) with two scales. In the left plot, the point shown in the filled circle at t = 0.039 represents the local minimum 
of ipjfc) and leads (p, a) = (0.87626, 0). This reduces to the one-parameter case. The filled circle at t = 0.146 gives negative a 
and we neglect the solution. In the right plot, we emphasized two candidates both of which provide positive two parameters. 
The gray box indicates the element of good sequence of /3 C estimation. 



TABLE II: Estimation of /3 C , pi and pi with two parameters. 



order 20 21 22 23 24 25 

Pc 0.221442 0.22140 0.221505 0.221513 0.221553 0.221562 

1/pi 1.138882 1.127291 1.160289 1.163069 1.179823 1.183677 

1/pa 0.433345 0.413636 0.468705 0.472926 0.500057 0.505965 

t* 0.116906 0.096589 0.119071 0.114190 0.120643 0.116991 



the transit region from high temperature to scaling regions. The right region separated by the right vertical dotted line 
is also not physically interesting, because in the re gion, a becomes negative. Even if we would admit small negative 
a for the compensation of a bit too large p- value, (3.26) gives 



P = 



1 1 

— + — • 

Pi P2 
1 

PlP-l 



(3.35) 
(3.36) 



and negative a* leads to negative P2 for pi > 0, which contradicts to < pi < P2 < ■ ■ ■■ Hence, the case of negative 
a* does not provide us truly reliable estimation of p\ and pi- From these consideration, we conclude that the scaling 
region is roughly implied by the two dotted vertical lines and only the cases a* > are worth of serious consideration. 
In the right plot in Figure 5, two such solutions labelled by black triangle and gray box are shown. 





1 


t 1 















Pc(25) 



0.05 0.1 . 15 0.2 0.25 



FIG. 6: ip2s{p,or,t) for p = 1.6896423, a = 0.5988992. The two dotted lines indicate /3 C = 0.22165 and its estimation, 
&(25) = 0.221562. 

The result of estimation is shown in Table 2, Figure 6 and Figure 7. The behavior of -ijj^ becomes steady at 12th 
order and single solution is obtained up to 19th order. The new channel to the most accurate sequence opens from 
20th for even orders and from 23rd for odd orders. This sequence includes estimation lying in the inside of scaling 
region. The level of the scaling behavior observed in ^25 is, however weak compared to the one-parameter case. The 
scaling level of 25th order in two-parameter case is, to the eye, the same level with the 13th order in one-parameter 
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0.22225 i 

0.222 
. 22175 Pc 

. 2215 
. 22125 

12 14 



16 18 20 22 24 

order N 



2v 



18 20 

order N 



FIG. 7: Plots of estimated f3 c and p 1 1 = 2«/ with two-parameters. The dotted lines indicate /3 C = 0.22165 and 2v = 1.26. The 
gray boxes indicate accurate sequence of which elements appears from 20th order. 



case. However the accuracy is improved. The relative error of 25th order (3 C is about 0.004 %. On the other hand, 
the estimation of p^ 1 — 2i> is not so good, though two-parameter estimation is improved compared to one-parameter 
case. At 25th order, the relative error is about 6 %. 

The extrapolation of the good sequence to the N — > oo limit is not adequate in two-parameter case. This is because 
the number of elements in the sequence is not enough. 



D. /3 C and v in three-parameters case and more 

With three parameters p, a and r, we deal with 

*P(p, a, r; x) = [l + pxj- + o\x^-f + r(xj-f] f3. (3.37) 

The estimation procedure follows those of one- and two-parameter cases and we omit the details and present just 
the outline. The stationarity condition reads ipff = 0, (i = 1,2,3) and then p, a and r are given by the following 
equations as function of t, 

^ + # 1) +# 2) +# 3) =0, = 1,2,3). (3.38) 

Then, consider the minimization of anc ^ the determination oft* gives optimal set (p*,cr*,r*) and 

p c (N)=$ N (p*,a*,T*,t*). (3.39) 
Estimation of pi (i = 1, 2, 3) can be given by (p* , a* , r*) via 

111 

9 = — + — + —, 

Pi P2 P3 

111 

a = h 1 , 

P1P2 P2P3 PzPi 

t = — 1 — . (3.40) 
P1P2P3 

With three parameters, ?/>jy shows complicated behavior to 22th order. Only from 23th order, our method begins 

to provide estimations characteristic to the three-parameter case. Even then, the behavior of ipff is not matured yet 
compared to the higher orders in the one- and two-parameter cases. It is typically reflected to the narrowness of the 
scaling region where all of p, a and r are positive. In addition, at the last even order TV = 24 in three-parameter 
case, the behavior of if: 1 ^ resembles to the 12th or 14th orders in the two-parameter case. See the left plot in Figure 
8. Table 3 shows the results at 23rd, 24th and 25th orders. As results in the first several orders exceeded of the 
correct value of /3 C in one- and two-parameter cases, these three estimations exceed (3 C — 0.22165, though the last 
order estimation is most close to the established value. We notice that at 25th order, the estimation of p^ 1 is also 
most accurate. This implies that when j3 c is precisely obtained, estimated v is also accurate. The opening of the 
accurate sequence for the three-parameter case demands further computation of high temperature expansion, maybe 
up to 30th order or more. 
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TABLE III: Estimation of /3 C and pi and p2 with three parameters. 



order 


23 24 


25 




0.222079 0.221741 


0.221687 


1/Pl 


2.010008 1.336976 


1.274452 


1/P2 


0.920863 0.715954 0.648486 


1/P3 


0.299087 0.242590 


0.210993 


t* 


0.111805 0.111116 


0.112973 



Pc(25) 
Pc 



¥ 25 (p,a,T,t) 

0.05 0.1 0.15 0.2 0.25 0.05 0.1 0.15 0.2 

FIG. 8: Left plot: (t). There are five points at which Wl becomes locally minimum. The left gray circle gives positive 
p, a and r. However, the two of pi (i = 1,2,3) becomes complex conjugates. The two black circles give r = and r < 
solution, respectively from the left to the right. The solution r = corresponds to the solution in the 2-parameter case with 
smaller t — 0.0837 (Represented by a black triangle in the right plot in Figure 5. This solution lies outside of scaling region and 
not included in the good sequence.) The rightest candidate at t = 0.2067 gives two of p; (i = 1,2,3) negative and we neglect 
it. Thus, the possible choice is only the black triangle for which f3 c and pi give values shown in Table 3. The region where the 
black triangle included is too narrow and this implies that the order of expansion is still short for clear demonstration of the 
scaling for 3-parameter ip N . Right plot: ii 2 s{p,cr,T,t) for p = 2.1339308, a = 1.2321906 and r = 0.1743783. The dotted line 
indicates /3 C (25) to the established digits, j3 c = 0.22165. 
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Under the simple estimation without using extrapolation to the infinite order, increasing the number of free pa- 
rameters improves the estimation of critical quantities so far. However, the reliable estimation with confidence of the 
scaling behavior of relevant functions sets in larger orders when the number of parameters are increased. This stems 
from the fact that the differentiation on f3^(t) creates oscillation and delays the appearance of scaling behavior. To 

say in the detailes the rationale is as follows: The small t behavior of P^\t) reads 

oo 

ffi® ~-J2 A " C ( N >-Pn)(-Pn) k t- P "- (3-41) 
n=l 

From < pi < p2 < P3 < ■ ■ ■ and the result in Table 3, we find p 2 > 1. Then for n > 2, p k n (n > 2) grows with k. This 

-(h) 

means that the differentiation enhances higher order corrections and the critical behavior of fi N is obscured. While 
at small t, the differentiation on /3 creates n to the coefficient of t n . Hence, also in small-i expansion, the upper limit 

~(k) 

of effective region of j3 N tends to shrink. We have actually found that, up to 25th order, m-parameter extension 
when m > 4 does not work. In our study up to 25th order, the three-parameter case is at the limit of multi-parameter 
extension of j3 being effective. 



IV. DISCUSSION AND CONCLUSION 



It would be better to mention on the P2 estimation. As would be understood from Table 2 and Table 3, the 
estimation of p2 is not successful in our approach. If one uses standard values for v = 0.63 and 9 — 0.5, one has 
p^ 1 ~ 0.84. Our estimated results in two- and three- parameters are still far from the value. Perhaps, P2 and also pi, 
may be obtained with more accuracy if we can invent the method where we directly address to pi estimation. In our 
approach, parameters p, a and r are just optimized in order to estimate /3 C . 
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Now, let us give another comment. The behaviors of the sequences of (3 C and v share among themselves a similar 
pattern according to the expansion order N. The pattern is such that the estimated result decreases first and, after 
the bottom out, turn to increase and then the good sequence appears (see Figure 4 and Figure 7). Also when the 
number of parameters was increased from single to double, same pattern repeated though the onset of the pattern 
delayed for several orders. The three parameter case would also follow this course. On the sequence of j3 c (N) first 
appearing, we note that any element stays within pre-scaling region. It might shoot wrong value in the TV — > oo 
limit. In fact, the extrapolation in the simple ansatz used in section [ill B| yields value plagued with non-negligible 
discrepancy with the correct value. Presumably the three parameter case is not exceptional, too. 

With focus on the scaling behavior of ip ~ Sfe=o Pkfi^ k \ we have attempted the estimation of /3 C and v. In the 
one- and two-parameter cases, good sequence emerged in the timing that the scaling regions of all relevant functions 
including ip, ' ' ' begin to appear. That sequence is most important since the element lies in the deepest place of 
the observed scaling region. The accuracy of the estimated results of our approach is not so superior compared with 
other approaches explained in [7]. Especially, the estimation of v needs more improvement for the higher accuracy. 
This remains as a problem of our approach. However, in our approach, it is transparent that which one is important 
when two or more candidates are present at an order. It becomes possible because the critical behaviors in relevant 
functions can be directly visible. At a fixed order, we can thus find single best estimation systematically, which is 
desirable for theoretical approach. 
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